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ABSTRACT 


By  the  calculations  presented,  the  minimum  total  velocity  increment  required 
for  bi -elliptic  transfer  between  non-coplanar  circular  orbits  is  obtained.  The 
maneuver  considered  is  the  following:  A  vehicle  in  circular  orb'.t  at  altitude 
h.  (radius  r.)  applies  an  impulsive  velocity  AVj  at  the  line  of  nodes.  The 
effect  of  the  application  of  AV^  causes  a  plane  change  of  amount  and  a 
transfer  ellipse  to  a  given  transfer  altitude  h^  (radius  r^)  is  established. 

When  the  vehicle  reaches  h^,  a  second  impulsive  velocity  change  AV^ 
simultaneously  changes  the  plane  by  amount  and  initiates  a  transfer 
ellipse  to  the  altitude  h^  (radius  r^)  of  the  target  orbit.  A  last  impulse 
changes  the  plane  by  amount  and  circularizes  the  orbit  at  altitude  h^, 
placing  the  vehicle  in  the  final  (target)  circular  orbit. 

Studies  were  made  of  the  choice  of  plane  change  angles  a^,  a^,  and  a^,  which 
minimizes  -  AVj  +  A  +  AV,  for  given  values  of  In,  h^  h^  and  total 

plane  change  angle  0  =  aj+a2  +  a^-  Several  limiting  relations  were  obtained 
for  Qj,  Q2«  and  a^I  they  are  dependent  on  either  the  ratio  r^/r.  alone,  or  the 
ratios  rt/r.  and  rf>  and  are  independent  of  0. 
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I.  INTRODUCTION 


Problems  dealing  with  orbital  transfer  are  of  considerable  current  importance. 
For  non-coplanar  orbits,  it  is  particularly  important  to  minimize  (or  nearly 
minimize)  the  velocity  expenditure  necessary  to  accomplish  the  specified 
plane  change.  In  this  report, the  minimum  velocity  expenditure  for  bi -elliptic 
transfer  between  non-coplanar  circular  orbits  will  be  obtained. 

One  use  of  bi -elliptic  transfers  is  in  the  physical  problem  where  rendezvous 
is  desired  between  vehicles  in  non-coplanar  circular  orbits.  The  line  of  inter¬ 
section  of  the  two  orbit  planes  will  be  referred  to  as  the  line  of  nodes.  The 
altitudes  above  the  surface  of  the  earth  of  the  inner  and  outer  orbits  are  re¬ 
spectively  denoted  as  In  and  h^.  and  the  plane  change  angle  is  6. 

In  this  study,  plane  change  angles  in  the  range  0  <  0  £  ir/2  will  be  considered. 
It  should,  however,  be  noted  that  the  various  conclusions  relative  to  the  divi¬ 
sion  of  plane  change  angle  and  the  equations  for  obtaining  plane  change  angle 
remain  valid  for  ir/2  <  0  £  it.  The  interest  in  plane  change  angles  greater 
than  tt/2  did  not  appear  sufficient  to  warrant  the  necessary  additional  length 
and  complexity  in  the  various  proofs.  A  study  of  bi -elliptic  transfers  with 
0=0  can  be  found  in  Reference  1. 

The  bi-elliptic  transfer  is  initiated  by  applying,  at  the  line  of  nicies,  a  velocity 
increment  AVj,  which  transfers  the  rendezvous  vehicle  from  the  circular 
orbit  at  altitude  In  into  an  elliptical  orbit  with  apsidal  altitudes  In  and  h^, 
while  simultaneously  rotating  the  orbit  plane  through  an  angle  a^. 

A  second  velocity  impulse  AV^,  applied  at  altitude  h^,  simultaneously  trans¬ 
fers  the  rendezvous  vehicle  into  a  second  elliptical  orbit  with  apsidal  altitudes 
h^  and  h^,  while  rotating  th-?  plane  through  angle  a^,.  The  rendezvous  maneuve 
is  completed  by  the  application  of  velocity  increment  AV^,  which  transfers 
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the  rendezvous  vehicle  from  -the  second  elliptical  orbit  into  the  final  circular 
orbit.  Evidently  the  final  plane  change  angle  is  =  0  -  -  a^.  Figure  1 

is  a  sketch  of  the  bi -elliptic  transfer  maneuver  when  r^  <  r^  <  r^. 

For  a  rendezvous  misaion,  the  choice  of  h^  depends  upon  the  relative  phasing 
of  the  two  orbital  vehicles.  References  2  and  3  discuss  the  selection  of  h^ 
for  a  given  rendezvous  mission  and,  in  addition,  present  an  alternative  three- 
dimensional  transfer  scheme.  Since  ht  is  determined  by  methods  described 
in  the  references,  the  present  analysis  can  be  considered  a  treatment  of  a 
pure  transfer  problem  in  which  h^  is  given  along  with  h.,  h^,  and  0. 

sf 

Three  separate  cases  must  be  considered,  which  can  be  classified  as  follows: 

h.  <  hj  ^  ht  (outer  bi -elliptic  transfer) 

h.  S  h  <  h  (intermediate  bi -elliptic  transfer) 

1  »  I 

h^  <  In  <  h^  (inner  bi -elliptic  transfer).  (1) 

The  velocity  required  to  transfer  from  h^  to  h^  to  In  is  equal  to  the  velocity 
necessary  to  make  the  transfer  in  the  reverse  order.  Therefore,  the  three 
cases  in  Eq.  (1)  cover  all  logically  possible  cases,  with  the  possible  relabel¬ 
ing  of  In  and  h^.. 

It  is  desired  to  determini  what  portion  of  the  total  plane  change  0  should  be 
accomplished  at  each  of  the  altitudes  concerned  (In  ,  ht>  and  h^)  in  order  to 
minimize  the  total  velocity  expenditure.  The  plane  change  angles  at  In,  ht> 
and  hj  are  respectively  denoted  by  ,  a^,  and  a^.  The  corresponding 


The  degenerate  case  In  =  h^  is  not  treated  here. 
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Ellipses 


Figure  1 


The  General  Bi-elliptic  Transfer  Maneuve 


geocentric  radii  are  defined  as 


r.  =  r  +  h.  ,  r  =  r  +  h.  ,  r-  =  r  +  h.  (2) 

iOl  tot  fof 

* 

where  rQ  is  the  radius  of  the  attracting  sphere. 

Section  II  describes  the  mathematical  technique  to  be  utilized  in  the  subsequent 
analysis.  Pertinent  results  are  presented  and  briefly  discussed;  additional 
graphical  results  are  presented  in  Section  VII. 


For  present  purposes,  the  attracting  sphere  will  be  assumed  to  be  the  earth. 
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II.  SUMMARY  OF  RESULTS 


Bi-elliptic  transfers  with  plane  changes  whose  total  velocity  requirements 
are  minimized  all  possess  two  (of  the  three)  angles  that  are  bounded.  For 
the  outer  bi-elliptic  transfer,  they  are  the  first  and  third  plane  change  angles 
( a j  and  a^).  For  the  inner  bi-elliptic  transfer,  they  are  the  first  and  second 
angles  (a^  and  a^).  Either  and  a,,  or  and  a^,  are  boarded  for  the 
intermediate  bi-elliptic  transfer. 

A  numerical  upper  bound  can,  in  fact,  be  placed  on  some  of  the  above  angles 
as  follows: 

a.  The  angle  a}  is  less  than  an  angle  K,  which  is  approxi¬ 
mately  5.  30  deg  for  either  the  outer  or  the  intermediate 
bi-elliptic  transfer. 

b.  For  an  outer  bi-elliptic  transfer  (13  <  K  =  5.  3  deg. 

c.  For  an  inner  bi-elliptic  transfer  <  K  -  5.  3  deg. 

The  other  angles  cited  are  bounded  by  functions  of  the  orbital  radii  and  the 
transfer  radius. 

The  Hohmann  transfer  is  inT  -  agated  as  a  limiting  case  of  a  bi-elliptic 
transfer. 

If  Qj,  a^,  and  are  respectively  the  plane  change  angles  at  r.,  rt,  and  r^, 
then  the  velocity  increment  AV,p,  necessary  for  the  bi-elliptic  transfer,  can 
be  expressed  as 


AV 


T 


V~ 

ci 


f(x,y,a1,a2,a3) 


where  V  .  is  the  circular  orbit  velocity  at  radius 


and 


x 


f 


>  ri.) 
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The  total  plane  change  0  is  simply  the  sum  Qj  +0,2  + 


The  values  of  a^,  a^,  and  a^,  which  minimize  AV,p/Vc.,  are  respectively 
denoted  as  a.  ,  a_  ,  and  a,  .*  If  r  2  r4,  the  angles  a.  and  a-  are  ob- 

18  m  8  jS  t  I  10  J  6 

tained  as  the  unique  solutions  of  the  equation 


=  0 


(3) 


The  following  properties  of  and  can  be  deduced  from  Eq.  (4),  when 


a,  <  a,  <  K  =  5.  30* 
Is  1 


a-  5  a,  <  K 
38  3 

Qj  <  for  all  y  <,  x 

Qj  =  cij  for  y  =  x  (r^  -  *f) 

Since  a,  +  eu  +  a,  =0,  any  two  of  the  above  three  angles  determine  the 

,1  •  1  10  C  8  J  8 

third. 

All  angles  are  between  0*and  180*  or  0  and  it. 


-6- 


Qls  =  Q1  =  ^  ’  when  x  =  1 

Qj8  =  =  0,  when  y  =  1 


when 


x  =  1  ,  xiy  =  1 

and,  therefore, 

a^fl  =0,  when  x  =  1 

The  properties  when  x  =  y  =  1  can  evidently  be  deduced  on  purely  physical 
grounds.  The  case  x  >  1,  y  =  1  is  seen  to  be  the  Hohmann  transfer  with  plane 
change.  The  above  result  (a3g  =  0  when  y  =  1)  implies  that  no  savings 

can  be  achieved  by  making  part  of  the  total  plane  change  after  circularizing 
at  the  final  orbit  altitude. 

The  subsequent  study  of  the  case  x  =  1,  y  <  1  indicates  that  no  plane  change 
should  be  initiated  prior  to  the  entrance  into  the  transfer  ellipse.  Therefore, 
for  a  Hohmann  transfer  the  total  plane  change  maneuver  should  be  divided 
between  the  initial  impulse  removing  the  vehicle  from  the  inner  orbit  and  the 
final  impulse  placing  the  vehicle  in  the  outer  orbit.  Furthermore,  the  initial 
plane  change  angle,  resulting  in  the  minimum  value  of  for  a  Hohmann 

transfer,  is  less  than  a.j. 

When  rt  <  r.,  it  is  found  advantageous  to  solve  Eq.  (3)  for  and  a^, 

rather  than  a,  and  a,  .  If  r.  S  r.,  then  a.  <  aC  and  a,  £  al  where 
Is  3s  t  i*  Is  1  Zb  2 
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The  following  properties  of  and  can  be  deduced  from  Eq.  (5),  when 

rtsV 


a,  <  a'.  S  a,  =  arc  cos 
Is  1  lm 


x  f 

«1.  =ai  =°.  when  -  =  —  =oo 

y  i 

Qj  =  alm,  when  y  =  x(r.  =  rf) 

a.  =  a'.  =  0,  when  x  =  l(r.  =  r.) 

Is  1  *  t  l 

a,  Sa'  <K  -5.  30* 

2s  2 


t 

a.  =  ai  =  0,  when  x  =  —  =  0 

2s  2  *  r. 

l 

a2s  =  a2  =  °*  when  Y  =  x(ri  =  rf) 

a'  is  maximized,  when  x  =  l(r  =  r.) 

u  t  1 

It  is  interesting  to  note  that  the  angles  a?  ,  j,  and  a £  are  all  bounded  by  the 
same  angle  K  =  5.  30*. 

For  r.  <  r.  <  r,,  a.  s  ei  (see  Eq.  (4).  Furthermore,  when  r.  <  r.  <  r- 
ltf’lsl  ^  *  ltf 
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if 


1  +  y  , 
x  >  — -f-  -  1 

y 


Similarly, 


Q3s  <  a3m 


if 


y 

The  present  analysis  places  bounds  on  the  independent  variables  regardless 

of  the  values  of  r^,  rt>  and  r^.  The  two  variables  chosen  (from  the  three 

variables  a.  ,  ,  and  a-  )  are  in  turn  determined  by  the  values  of  x  and  y. 

18  68  ^8 
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HI.  TRANSFER  ALTITUDE  GREATER 
THAN  FINAL  ALTITUDE  hf 


Consideration  will  begin  with  the  first  classification  in  Eq.  (1).  with  which, 
in  summary,  the  following  mane  we  r  is  associated.  A  transfer  is  initially 
made  from  a  circular  orbit  of  altitude  h^  to  ail  elliptic  orbit  with  apogee 
altitude  ht>  The  angle  between  the  two  orbit  planes  is  a^. 

It  can  be  shown  that  the  velocity  increment  AVj,  necessary  to  accomplish 
this  initial  maneuver,  is  given  by 


AV1  I - 1 

-V~L  =  Jl  -2HlCos  Qj  +  Hj 

ci 


(6) 


The  functions  H.  and  V  .  are  defined  as 
1  ci 


2x 
+  x 


(7) 


where  x  =  r./r.,  u  is  the  force  constant  for  the  earth,  and  the  function  V  . 

t  1  r  ci 

is  the  circular  orbit  velocity  at  altitude  h.. 

The  second  transfer  maneuver  is  initiated  at  altitude  ht>  It  consists  of  trans¬ 
ferring  from  an  orbit  with  apogee  at  ht  and  perigee  at  h^  to  one  with  apogee 
at  ht  and  per. gee  at  h^,  accompanied  by  a  plane  change  at  h^. 

The  velocity  A  corresponding  to  the  second  maneuver  is  given  by 


AV 


2 


V  . 

ci 


-  2H^cos  a ^ 


(8) 
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The  functions  and  are  defined  as 


H2  =JxTT  +  yT 


(9) 


where  x  =  rj r.  and  y  =  r^/ r^. 

The  final  maneuver  consists  of  transferring  from  an  elliptical  orbit  with 
apogee  at  and  perigee  at  h^  to  a  circular  orbit  at  altitude  h^,  accompanied 
by  a  plane  change  a^.  The  velocity  expenditure  AV^  can  be  given  as 


where 

=  ^  •  H5  *  anda,  =e-0l-a2  (It) 

It  follows  from  the  above  analysis  that  the  velocity  increment  AV^,  for  the 
bi -elliptic  maneuver  is  given  by 


where  the  three  components  are  given  in  Eqs.  (6)  through  (11).  It  is  desired 
to  obtain  the  minimum  AV,j,  for  given  values  of  x,  y,  and  0.  Equations  (6) 
through  (12)  can  be  shown  to  apply  for  arbitrary  values  of  h^,  h^,  and 
ht(h.  <  hf). 

It  will  be  shown  that  the  minimum  value  of  AV,p  can  be  obtained  from  an 
examination  of  the  partial  derivatives  of  AV,p  with  respect  to  and  a^.  A 
method  will  be  developed  to  obtain  this  minimum. 
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c 


From  Eqs.  (6),  (8),  (10),  and  (12),  the  partial  derivatives  of  with 

respect  to  and  a^,  are  given  by 


1  8avt 

Vci  "S 


a 


m 


HjSin  Oj 


H2H3*in  °2 


il  -  2HjCos  Qj  +  Hj  Jl  -  2h7co*  a  +  Hj 


1  3*VT 


H,H_sin  a, 
4  5  3 


H2H38in  a2 


§a, 

ci  3 


I 


1  -  2H^cos 


ZH  cos  +  Hj 


(13) 


where  =  0  -  -  a^. 

Examination  of  Eq.  (13)  shows  that  (when  y  s*  x  and  y*  1): 


9AVt 

-B5T  <0at  °l  =° 


3AVt 

>  0  at  at  =  0  -  q3  ; 


9AVt 

”55^-  <  0  at  q3  =  0  , 


aAVT 

“SoT  >  0  at  a3  =  6  •  ai 


* 


(14) 
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The  first  pair  of  values  in  Eq.  (14)  implies  that,  for  any  a^(0  i  ^  6),  there 
exists  an  a^,  which  minimizes  AV^,  for  the  chosen  value  of  a^.  It  further 
follows  that 

9  A  VT 

f(al*a3)  = — Sa^~  =  0  (15) 

at  the  above  minimum;  i.  e.  ,  the  solution  of  Eq.  (15)  yields  the  value  of 
which  minimizes  AV^,  for  any  chosen  value  of  a^.  Similarly,  the  second 
pair  of  values  in  Eq.  (14)  implies  that  there  exists  an  a^,  which  minimizes 
AV^,  for  any  given  a^(0  S  <  0).  The  locus  of  these  values  is  obtained 
by  finding  the  roots  of  gfa^cij)  =  dAV^/da^,  i.  e.  ,  by  solving  the  equation 

9AV„ 

«<Vj»— =  °  <16» 

for  any  given  value  of  a^.  The  minimum  value  of  AV^,  is  obtained  by  solving 
the  following  system  of  equations  for  and  (see  Eqs.  (15)  and  (16)): 

*(«*!.  a3)  =  0 

gtej.Qj)  =  0  (17) 


if  a  solution  exists. 

Physical  reasoning  indicates  that,  for  h^  >  h^,  the  values  of  Qj  and  which 
minimize  A  V,p,  should  be  small.  This  fact  alone  is  usually  sufficient  fo."  ob¬ 
taining  iterative  solutions  to  Eq.  (17)  (see  Reference  4).  The  investigation 
that  follows  provides  further  insight  into  the  problem  of  minimizing  AV^.  for 
all  values  of  h.,  h  ,  and  h-(h.  *  h,). 

I  t  I  X  I 
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The  exigence,  uniquene...  and  aingularitie.  of  .olution.  of  £,.  ,17)  ^11  ^ 

d;acu  d  .n  Sec, ion  fV.  fn  order  to  accotnpfiah  ,,  lnve.h^io 

attenhon  ,s  f.r.t  directed  to  the  function  V' 


F(a)  = 


H  sin  a 


/l  +  H2  - 


2H  cos  a 


(18) 
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IV.  PROPERTIES  OF  THE  FUNCTION  F(a) 


The  following  properties  of  F(a)  are  deduced  from  inspection: 


F(a)  >  0  for  H  2  0 


F(0)  =  0  for  H  *  1 


F(i)  = 

VI  +  H 


a  9) 


Setting  the  first  derivative  of  F/n)  equal  to  zero  shows  that  the  maximum 

value  of  F(q)  occurs  when  a  =  a  where 

m 


1  +  H  . 

co* 


2  2 
/  1  +  H  r 

\l\  ZH  /  ~ 


1  +  H 

2H 


y/d  -  H2)* 
2H 


(20) 


Since  ( 1  +  H2)  / 2H  2  1  for  all  H,  it  follows  that  the  maximum  value  of  F(a) 


occurs  at  a  =  a  where 
m 


cos  a  =  'or  K  21 
m  rl 


cos  a  =  H  Ijt  H  <  1 
m 


(21) 


Substitution  of  Eq.  (21)  into  Eq.  (18)  yields  the  maximum  vaue  of  F(n)  as 

F  =  1  for  H  >  1 
max 


F  =  H  for  H  <  1 
max 


(22) 
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V.  SOLUTION  OF  THE  EQUATIONS  FOR  THE 
MINIMUM  TOTAL  VELOCITY 


This  section  will  reformulate  the  system  of  equations  in  Eq.  (17)  (Eqs.  (15) 

and  (1 6))  in  a  form  more  suitable  for  iterative  solution.  The  task  of  iteratively 

solving  Eq.  (17)  will  be  further  facilitated  by  obtaining  the  upper  bounds  on  the 

angles  a.  =  a.  and  a0  =  a.  ,  which  minimize 

lls  ojs  i 

Substitution  of  the  first  equation  of  (13)  into  Eq.  (15),  when  8AV^,/9aj  =  0,  yields 

Hj.in  dj  H2H3sina2 

Gfap  =  - . .  . . =  — -  - - 

J 1  -  ZHjeos  a  +  Jl  -  2H3cos  a2  +  H^ 

H_H_sin(0  -  a.  -  a~) 

=  H(q2)  =— .  -----  —  =  F(a1>a3)  (23) 

Jl  -  2HjCos(0  -  Oj  -  Qj)  + 

From  Eqs.  (7)  and  (9) 

1  <  Hj  S  2 
0  <  H3  -  1 

0  <  H2  <  1  (24) 

It  follows  from  Eqs.  (22)  and  (24)  that  the  maximum  value  of 

H.  sin  a. 

g(q1>  =  ---  . .  1  -  -■ 

y  1  -  2H^cos  +  Hj 
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is  unity.  Similarly,  the  maximum  value  of 


H(a2)  = 


H2H3*in  °2 


l ~ 


2HjCoe  a2  + 


is 


H2H3 


=  v/Wrr 


5T 


<  1  (see  Eq.  9) 


From  Eqs.  (18)  through  (21)  (and  the  above  discussion),  the  function  G(a.j) 
increases  monotonically  from  zero,  when  =  0,  to  unity,  when  =  Qjm» 

It  then  decreases  monotonically  with  for  >Qim*  The  function 
H(a2)  =  H(6  -  Qj  -  a^)  increases  monotonically  with  from  H(0  -  a^)  £  0, 
when  cij  =0,  to  its  maximum  value  ls/T7x(l  +  x)  <1,  when 

a,  =  a  (a.  =  0  -  a,  -  a_  ).  The  function  I  ,)  then  decreases  monotonically 

C  cm  i  j  Cm  c 

with  from  the  maximum  value,  when  =  0  -  -  a2m»  to  zero>  when 

Qj  =  0  -  a3(a2  =  0)*  The  variation  with  of  G(Qj)  and  Hfa^)  is  shown  in 

Figure  2.  It  follows  from  the  above  that  Eq.  (23)  (the  first  in  equation  (17)) 
must  have  a  solution  for  a.j  <  (see  Eq.  (21)),  where 


a.  =  arc  cos 
1m 


arc  cos 


I 


1  +  x 

~ET 


(25) 


as  is  illustrated  in  Figure  2.  The  possibility  of  additional  solutions  when 
Oj  >aim  is  investigated  below. 

As  novid  above,  G(Qj)  decreases  monotonically  with  for  >  Qjm»  reach¬ 
ing  its  minimum  when  =  0  -  a^.  It  therefore  follows  that  Eq.  (23)  can  have 
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no  solution  for  a,  >  a,  if, 
1  1  m 


min  G(cij)  = 
al  >  Qlm 


(see  Figure  2) 

H.  sin(0  -  a.) 

1  J  > 

J 1  -  ZHjCoste  -  a3)  +  Hj 


H-H-sin  a_ 

maxH(a^)  =  max - 

J 1  -  ZH^cos  Qj,  +  Hj 


*  H2H3  * 

where  is  a  function  of  x  given  by  Eq.  (7).  Equation  (26)  will  be  shown  to 
be  true  over  all  but  a  narrow  band  of  possible  x,  0-values. 

Equation  (26)  is  true  if,  and  or.ly  if 


=  arc  cos  x 


(27) 


In  summary,  it  is  impossible  for  Eq.  (23)  to  have  a  solution  for  >  ajm»  if 
Eq.  (27)  is  satisfied.  Since  the  present  analysis  is  restricted  to  plane  change 
angles  in  the  range  0  S  90 *,  there  can  be  no  solutions  of  Eq.  (23)  for  Qj  > 
if 


cos  6  <0  (28) 

m 

Since  Eq.  (27)  must  be  true  if  Eq.  (28)  applies. 
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From  Eqs.  (27)  and  (28),  it  follows  that  for  any  0  S  90*,  Eq.  (23)  cannot  have 

a  solution  for  a,  >  a,  if 
1  lm 


f(x)  =  x^  + 


3x  -  1  >  0 


(29) 


Equation  (29)  is  true  if,  and  only  if 

x  >  x  =  1.48 


(30) 


where  x  is  the  single  positive  root  of  the  cubic  equation  f(x)  =  0  (see  ^.q.  29). 

Furthermore,  x  =  r./r.  2  r,/r.. 

t  l  f  l 

Figure  3  is  a  graph  of  0^  versus  x.  It  will  be  assumed,  without  proof,  that 
any  solution  g  Eq.  (23)  which  minimizes  AV^.  for  a  given  value  of 

satisfies  the  inequality  in  the  narrow  x,  0  -  region  excluded 

from  the  above  proof  (that  area  above  the  0  versus  x  curve  in  Figure  3). 
Therefore,  the  values  of  tij  =  a^,  which  minimize  AVT  for  given  values  of 
x,  0,  and  a^,  are  obtained  by  solving  Eq.  (23).  The  above  <ijg  values  that 
satisfy  Eq.  (23)  also  satisfy  the  inequality  aj8  -  aim* 

Equation  (23)  can  be  restated  as 


cos  a 


1  =  I^[s  ±  -  s(l  +  h9+ 


(31) 


2  2 

where  S  =  [F(aj,a^)]  =  [Hfa^)]  .  The  two  choices  of  sign  in  Eq.  (31)  corre¬ 
spond  to  values  of  greater  or  less  than  Since  the  present  investiga¬ 
tion  is  restricted  to  values  of  less  than  or  equal  to  only  the  positive 

sign  need  be  considered.  It,  therefore,  follows  that 


cos  a 


i  =i^[s  +  isi-s(1  +  Hi)  +  Ht] 


(32) 

(cont. ) 


-23- 


I 
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Figure  3.  Value*  of  x  and  0  Where  Proof  That  tij  5  a  Applie* 


iS-’jgEKJ 


sin  =  j|-  F<Qj  ,  Qj)  Jl  -  2HjCOS  +  Hj 


v  ***** 


From  Cqs.  (18)  through  (22),  a.  is  a  monotonically  increasing  function  of 
F(aj,a3)  =  H(a2)  =  G(a^)  in  the  interval  0  S  -  °im*  Furthermore,  all  solu¬ 
tions  <ij  =  Qjg  of  Eq.  (23)  lie  in  the  above  interval.  If  is  that  value  of 
at  which 


F(a1.a3)  =  maxH(a2)  = 


it  follows  from  the  above  monoticity  arguments  that  aig  -  for  x»  Y*  and 6. 


Since  when  F(a^  ,  a^)  =  H(a2)  =  H^H^,  it  follows  from  Eqs.  (7),  (9), 

and  (32)  that 


—  I  2  ,  (x  -  1)  /x  +  2 

1  Jx3(  1  +  X)  *  J  x 


sin  Qj  =  +  J 1  -  cos  Qj 


Figure  4  is  a  graph  of  versus  x.  It  is  interesting  to  note 


a,  <  a,  <  K  =5.  30° 
Is  1 


regardless  of  the  choice  of  x,  y,  and  0.  It  can  be  shown  that  is  maximized 
when  x  =  (1  +  n/7)/2  -  1.  82. 


The  second  equation  of  (14)  (Eq.  (13))  is  treated  below.  A  bound  is  similarly 
obtained  on  the  angle  =  q38  which  minimizes  for  a  given  value  of  a^. 


Substitution  of  the  second  equation  of  (10)  into  Eq.  (13)  yields  the  algebraic  equatio: 


H  .H,  sin  a. 
4  5.  3 


H2H3sm  a2 


=  F(a1,a3) 


j  1  -  2H5cos  q 3  +  Hg  Jl  -  2H3cos  a2  +  H3 
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Figure  4.  Bounding  Curve  for  Possible  Solutions  of  Equation  (23) 


iS  »'-  '*  • 

U.mg  the  same  line  of  reasoning  as  in  the  investigation  of  Eq.  (23),  «  ], 
readily  shown  that,  for  any  o  there  exists  an  o  <  o  .  which  satisfies 

Rn  Met  ...1 _  J 


Eq.  (35)  where 


=  arc  cos 


*  *~j t i*  •  4.~  f 

-  .•  •  jV ,  ;  .  'L  ■  w  ufnw  $  A 

--I*  ... 

Hr 1  arc  co»  Jtt*  •  • 


It  can  also  be  shown  that  Eq.  (35)  cannot  have  a  solution  for  o,  >  a  if 

3  3m 


H4H5 


1  +  H 


W*3*r  >h2h3  -  U rig- 


Y  .  1 

1  +  3y  1  +  x 

The  inequality  in  Eq.  (34)  is  true  if,  and  only  if 


y  >  y 


where  y  is  the  single  positive  root  of  the  cubic  equation 

(0y3  +  y2  -  37  -  1  =  0  (39) 

Figure  5  is  a  graph  of  y  versus  x/y  =  r^/r.. 

1.  foiiow.  from  Eq.  (39)  that  7  <  1  for  rf/r.  >  3.  Since  y  =  r./r.  II,  Eq.  ,37) 

“  alWayS  trUe  When  x/V  -  Vi  »  ’•  Furthermore,  it  is  shown  above  that 
Eq.  (35)  cannot  have  a  solution  for  Oj  >  o3m  when  Eq.  (37)  i,  true.  The  remain- 

mg  investigation  is  therefore  concerned  with  values  of  r./r.  in  the  range 
1  s  rf/r.  <3.  t  i  b 
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iilliil 


Figure  5.  Variation  of  y  with  r.f 


■'jnm 


For  y  S  y  (see  Eq.  (38)),  Eq.  (35)  cannot  have  a  solution  for  >  a^m  if 


H4H5sin(6  -  aj) 

fl  -  2H5cos(0  -  aj)  +  Hg 


>  max 


H2H3.m 


/l  -  ZH^cos  &2  +  Hj 


713  <«> 


The  inequality  in  Eq.  (40)  is  true  if  and  only  if  (see  Eqs.  (26)  and  (27)) 


0  -  a,  <  0'  =  arc  cosl 

1  m 


I  2(1  +  y)  _  ( 14  uy)  y  +  2(1  +  y)  -  (1 4  uy)(l  4  3y)y  | 

fy3(l  +  uy)2  J  y  (1  +  uy)  J 


where  u  =  x/y  =  r^/ r..  Figure  6  is  a  graph  of  0^  versus  r^r.  for  various 
values  of  y. 

It  will  be  assumed  that  for  all  values  of  x,  y,  and  6,  as  it  was 

previously  assumed  that  Qjs  <  aim  *or  ^  x  ^  ®- 

By  analogy  to  the  transformation  of  Eq.  (23)  into  Eq.  (32),  Eq.  (35)  can  be  writtei 


i 

co‘  “3  =  SJTL5  +  J®2  -  ^U1  +  Hs) +  HJHi 


5/  '  4  5 


sin  = 


F(a  j,  a3) 


1  -  2H5cos  03  +  Hg 


where  S  =  [F(oj,a3)J  . 
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From  Eqs.  (18)  through  (22),  a.  is  a  monotonically  increasing  functton  of 

J  “  '  i  ;<1 


FtQj.Qj)  =H(a2)  = 


HjHgSin 


fl  -  2H5cos  a3  +  Hg 


(see  Eq.  (35))  in  the  interval  0  <  £  Q3m*  Furthermore,  all  solutions 

q3  =  a3g  of  Eq.  (35)(or,  equivalently,  Eq.  (42))  lie  in  the  above  interval.  If  a^ii 


that  value  of  at  which 


F(a1(a3)  =  maxFtcij ,  a3)  = 


1*  3'  -  vx(l  +  x) 


then  it  follows  from  the  above  monoticity  property  that  S  for  all  x,  y, 
and  0. 


Since  when 


F(ai «  “3)  =  maxF(al  •  “3>  =  • 


it  follows  from  Eqs.  (9),  (11),  and  (42)  that 


cos  a.  =  /  -2*1  +  .y>  +  l1 1,  *y)  y  +  2(*  +  y?  -  +  *yH1  +  3i 

V  y  (1  +  uy )  >J  y  (1  +  uy) 


sin  =  /l  -  cos  , 


sre  u  =  x/y  =  r^/r..  Figure  7  is  a  graph  of  versus  y  for  various  values 


of  r^/r.. 
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Figure  7.  Bounding  Curve  for  Solutions  of  Equation  (35) 


It  is  interesting  to  note  (see  Eq.  (34))  that 


a3g  £  a3  <  K  =  5.30°  (44 

for  any  choice  of  x,  y,  and  9.  Comparison  of  Eqs.  (33)  and  (43)  shows  that 
cij  =  a3,  when  y  =  x.  Furthermore,  the  maximum  value  of  occurs  when 
x/y  =  r,/ri  =  1. 

From  Eqs.  (33)  and  (43) 


Q3a~°  asy~‘ 


a  -*  0  as  x  -»  1  (45) 

Is  '  ' 


Also,  since  x  >  y,  s  x  -*  1 ,  y  -»  1  and  thus 


a3g-*0asx-*-l  *  (46) 


> 
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VI.  TRANSFER  ALTITUDE  ht  BETWEEN  INITIAL 

ALTITUDE  h.  AND  FINAL  ALTITUDE  h, 

1  f 

The  analysis  of  the  bi-elliptic  transfer,  with  r^  >  r^,  treated  in  the  previous 

sections,  was  restricted  to  the  case  where  r.  *  r^  (except  for  the  degenerate 

case  r.  =  r^  =  r^,  where  Qjs  =  Q38  =  0).  Consideration  will  now  be  given  to 

the  case  where  r.  <  rt  <  r,. 

1  t  f 

The  velocity  AV^,  required  to  transfer  from  altitude  h.  to  h^  to  h^  with 

plane  change  6  in  given  by  Eqs.  (6)  through  (12),  where  AV^  is  the  velocity 

increment  at  h  =  h. ,  AV_  is  the  increment  at  h  =  h,.,  and  AV,  is  the  incre- 

i  2  t  3 

ment  at  h  =  h^.  It  thus  appears  that  there  is  no  difference  between  the  cases 

r.  <r,  <  r.  and  r.  £  r.  <  r,.  It  can,  in  fact,  be  shown  that  when  r.  s  r.  <  r-, 
lftitf  *  ’  ltf* 

the  minimum  value  of  AV^,  can  be  obtained  by  solving  Eq.  (17). 

However,  when  r^  <  rfc  <  r^,  y  <  1  and,  therefore,  the  various  bounding 

theorems  such  as  Eqs.  (33)  and  (43)  no  longer  apply.  It  is  also  convenient, 

when  employing  iterative  methods,  to  utilize  independent  variables  that  can 

be  bounded  in  the  manner  a  ^  and  a^were  bounded  in  the  first  problem(0  S  S 

0  S  <  a^m).  Bounds  will  therefore  be  sought  for  two  of  the  three  variables 

Qj,  a^,  and  q^,  and  the  existence  of  a  minimum  will  be  established.  Due  to 

the  length  of  the  necessary  exposition,  uniqueness  will  not  be  demonstrated. 

The  uniqueness  proof  is  of  essentially  the  same  nature  as  the  one  for  the  case 

r.  £  r,  £  r,  in  Reference  5. 
lit 

Suppose  that  the  two  independent  variables  are  chosen  as  and  a^.  It  then 
follows  that  the  two  partial  derivatives  of  AV,p  are  given  by  Eq.  (13).  If 
x  1  then  the  conclusions  in  Eq.  (14)  also  apply. 

From  Eqs.  (7),  (18),  (21),  and  (22),  the  function 

H.  sin  a. 

G(a  i )  =  —  =  -  (47) 

Jl  -  2H^cos  +  Hj 


-35- 


increases  from  zero  when  Qj  =  0,  to  unity  when 


=  a 


1  m 


=  arc  cos 


iS" 


(see  Eq.  (2 5))- 

Furthermore,  G(a.)  decreases  as  a,  increases  for  a,  >  a.,  .  From  Eqs. 

1  1  1  lm  n 

(9),  (18),  (21),  and  (22),  the  function 

H?H,sin  cu 

H(a2)  =  H(G  -  Qj  -  a3)  =  ——±-±  c  (48) 

Jl  -  2H3cos  q2  +  H3 


increases  from 


H2H3sin(0  -  q3) 

J 1  -  2H3cos(0  -  q3)  +  H3 


when  =  0,  to 

mixH(a2)  =  <  1  <49> 

when 

cos(6  -  Qj  -  a3)  =  =  co.  a2m  . 

The  function  Hfc^)  decreases  monotonically  with  Qj  from  its  maximum  value 
to  zero  when  aj  =  0  -  03. 


© 

From  Eqs.  (13),  (47),  (48),  and  (49),  there  exists  an  <  airn  a* 

which  8AV^,/8aj  =0. 

__  .  li  ' 

It  is  readily  shown  that  <  a,  (see  Eq.  (33))  and,  therefor'-,  all  solutions 
Q1  =  als  °*  He  ^e^ow  *ke  curve  aA  =  in  Figure  4. 

From  Eqs.  (11),  (18),  and  (22)  (also  see  Eq.  (10)) 

H^H,.sin  Qj 

maxK(a^)  =  max - 

Jl  -  2H5cos  a3  +  h| 


=  Y 


i 


2 

x(l  +  y) 


(50) 


It  can  be  shown,  employing  Eqs.  (47)  through  (50),  that  the  equation  9AV^,/8aj  = 
must  have  a  solution 


a.  =  a-  <  a-  =  arc  cos 
3  3s  3m 

(see  Eqs.  (13),  (23),  (49),  and  (50))  if 


y 


2 


Ti^y 


or 


*  1  +  y 

x-~rL 


(51) 


If  Eq.  (51)  is  true,  then  the  methods  of  proving  the  existence  of,  and  obtaining 
solutions  cij  =  Qjg  and  to  Eqs.  (15)  and  (16)  (yielding  the  minimum 

value  of  AV,p  for  given  x,  y,  and  0)  are  identical  to  the  methods  utilized  for 
the  case  h^  £  hr  <  h^  previously  analyzed. 
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(  I  • 


Regardless  of  whether  Eq.  (51)  applies,  Qj  £  -►  0  as  x  -*■  1.  It  can  be 

inferred  from  continuity  arguments  or  demonstrated  directly  that  assumes 

its  limiting  value  =  0,  when  x  =  1.  It  was  previously  shown  that  *  0 

when  y  =  1  (see  Eq.  (45)).  From  the  two  limiting  cases  (o^  =  0  when  x  *  1, 

=  0  when  y  =  1),  it  is  concluded  that,  for  a  Hohmann  transfer  with  plane 

change,  no  plane  change  should  be  made  prior  to  entering  or  after  leaving  the 

transfer  ellipse.  Furthermore,  since  the  Hohmann  transfer,  x  =  1  or  y  =  1, 

was  previously  treated,  the  present  investigation  can  be  limited  to  the  region 

h.  <  h  <  h,. 
i  t  f 

Because  convenient  bounds  cannot  always  be  placed  on  the  variable  a^,  con¬ 
sideration  will  be  given  to  the  independent  variable  pair  and  a^.  The 
partials  of  AV^,  with  respect  to  and  are  given  by  the  expressions  (see 
Eqs.  (6)  through  (12)): 


8AV.J,  H^sincij  H^H^sin 


From  Eqs.  (47)  through  (52),  it  follows  that  (see  Eq.  (51)),  when 

x>±4*-  1  (53) 

y 

♦here  exist  an  <aim  anc*  Q2  =  Q2s  <  Q2m*  respectively  satisfy 


=  0  for  any 


°2 


(54) 
(cont. ) 
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8AV, 


-g—  =  0  for  any  Oj 

w 


(54) 


In  summary,  it  is  seen  that  two  angles  are  bounded  regardless  of  the  values 
of  x  and  y.  If  Eq.  (51)  applies,  then  and  a3s  <  a3m‘  ^  ^1*  (53) 

is  true  (and,  therefore,  Eq.  (51)  is  false),  then  Ujg  <  a^and  a2s  <  a2m*  *n 
either  case,  a^g  £  Qj,  where  is  defined  in  Eq.  (33)  and  depicted  in 
Figure  4. 


It  can  be  shown  that  the  two  equations  in  (54)  have  a  single  solution  =  ajg* 


a,  =  a,  when  the  inequality  in  (53)  is  true.  Since  q_  =  0  -  a,  -a,  ,  and  sine 
Ct  Lt  S  J  8  18  CB 


the  two  equations  in  (54)  are  equivalent  to  the  two  equations  in (15)  and(16),  it 
follows  from  the  above  discussion  and  the  discussion  subsequent  to  Eq.  (51)tha 
Eq.  (54)  has  a  single  solution  a,  =  a,  ,  a,  =  a,  for  all  x  and  y  when  h.  <  li  <  h. 
It  was  previously  demonstrated  that  the  solution  obtained  yields  the  minimum 
value  of  AV^,  when  Eq.  (53)  is  false  (see  Eqs.  (14)  and(51).  It  is  demonstrated 
below  that  the  solution  also  yields  the  minimum  value  of  A  when  Eq.  (53)  is  ti 


From  Eq.  (52) 


3AV 

— 5 -  <  0  when  a,  =  0 

da,  1 


9AVt 

— s -  >  0  when  a,  =0  -  a, 

da.  1  2 


3AVt 

— 5 <  0  when  a,  =  0 

oa,  2 


3AVt 

—r -  >  0  when  a,  =  0  -  a, 

da,  2  1 


(55) 
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It  follows  from  Eq.  (55)  that  the  single  solution  of  Eq.  (54)  corresponds  to  the 

minimum  value  of  AV  for  h.  <  h.  <  h- 

T  1  t  f 

When  Eq.  (53)  is  true  (a^  s  ~  a2m^’  t^ien  Eqs.  (52)  and  (54)  can  be 

combined  to  yield  the  following  equations: 


where  S  =  [FUj.a^]2  =  [K(q3)]Z. 
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VII.  TRANSFER  ALTITUDE  ht  LESS  THAN 

INITIAL  ALTITUDE  h. 

1 


In  the  previous  sections,  the  bi-elliptic  transfer  has  been  analyzed  in  the  regioi 

h^.  <  ht  and  h.  s  h{  <  h^.  It  therefore  follows  from  Eq.  (1)  that  the  only  region 

to  be  investigated  is  h  <  h..  As  previously  noted,  the  velocity  increments 

^  1  I 

for  any  values  of  h.,  ht,  and  h^  are  given  by  Eqs.  (6)  through  (12);  die  two 

partial  derivatives  of  interest  are  given  by  (52).  The  conclusions  in  Eq.  (55) 

also  apply  in  the  present  region. 

In  Section  V  bounds  were  obtained  for  the  two  independent  variables  and 
a^.  Similar  bounds  will  be  obtained  in  this  section  for  the  variables  Oj  and  a. 

In  the  region  <  h^  <  h^  (see  Eqs.  (7),  (9),  and  (11)) 


Hj  <  1  ,  H3  <  1  .  H5  <  i  (57] 


From  Eqs.  (7)  through  (11),  (18)  through  (22),  (47)  through  (50),  and  (57),  it 
follows  that: 


a.  The  function  G(aj)  increases  monotonically  from  zero 

when  c.,  =  0,  to  G(o.)  =  H.  when  a.  =  a.  =  arc  cos  H.. 
1  ’  'll  1  lm  1 

For  Qj  >  G(Qj)  decreases  monotonically  as  Qj 

increases. 


For  a  given  value  of  a^,  the  function  K(aj)  increases  mom 
tonically  from  a  value  K'  >  0  when  Qj  =  0,  to  K(a^)  =H^H 


when  Qj  =  a^m  =  9  -  Uj 


=  arc  cos  H^.  The  function 


Kta^)  decreases  monotonically  from  its  maximum,  when 

a.  =  9  -  a_  -  a-  ,  to  zero  when  a.  =  0  -  a,. 

1  2  lm1  1  2 

For  a  given  value  of  ,  the  function  K(a^)  increases  from 

K(a.j)  =  K"  >  0  when  =  0,  to  K(a^)  =  H^H,.  when 


°2  *  6  -  °1  -  a3m- 
=  6  -  a,. 


It  then  decreases  to  zero  when 
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d.  The  function  H(a^)  increases  monotonic&ily  from  zero 
when  =0,  to  H(a2)  =  H^Hj  when  =  a2m  =  arc  co* 

For  a 2  >  a2m’  **(».,)  decreases  monotonically  as 
increases. 

From  (a)  through  (d)  and  Eqs.  (7),  (9).  and  (11),  it  follows  that  for 
ht  <  h.  <  hf(x  >  y) 

maxK(a3)  =  < 

maxG(a^)  =  1-lj  =  JrT^  <  maxH<°Z>  =  Jx(l\  x)  ,58) 

V  X 

From  Eqs.  (52)  and  (58),  it  follows  that,  for  any  a?,  there  exists  an 
cij  =  Qis  <  Qjm  at  which  (see  (a)  through  (d)  above) 

G(Ql)  =K(q3) 


or 


9AV_ 

up  ■ 0  <59> 

Similarly,  for  any  q^,  there  exists  an  =  a2g  <  at  which 

H(q2)  =K(q3) 


or 


(60) 
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From  Eq.  (55),  there  exists  an  <  ajm(°2*  <  aZn)'  Eq.  (591 

(Eq.  (60))  and  yields  a  relative  minimum  of  AV ^(aj.a.j).  &  cam  further  be  she 
that  only  one  value  of  dj  ~  a  ^  g  satisfies  Eq.  (59)  for  each  value  of  a  and,  sii 
larly,  for  each  value  of  a^,  there  is  only  one  =  a , a  that  satisfies  Eq.  (60). 

It  can  be  shown  (see  Reference  5)  that  the  simultaneous  solution  of  Eqs.  (59) 
and  (60)  yields  the  values  of  =  Qjg  and  ®2  =  a2s  w^'®h  minimise  AV^, 
where  als  <  Qjm  and  Q2s  -  a2m-  Equations  (59)  and  (60)  can  be  written  in 
the  more  explicit  form  shown  in  (56). 

From  Eq.  (58)  and  the  properties  (a)  through  (d)  preceding  it,  the  function 
H(a2)  (Eq.  (48))  increases  from  zero  when  a2  =  0,  to  maxHfa^)  when  =  a2  . 
Furthermore,  a2g  <  a2m*s  the  solution  of  Eq.  (60).  If  S  a2m  t*le  v%*ue  °* 
for  which  (see  Eq.  (60)), 

H(a2)  =  H(a£)  =  maxK(a3)  s  K  =  H4H5  (61) 

and  it  follows  from  Eq.  (58)  that  a2s  £  a2  for  all  x,  y,  and  0.  From  Eqs.  (9), 
(11),  and  the  third  and  fourth  equations  of  (56),  a2  is  given  by 


where  u  =  x/y  =  r^/r^  Figure  8  is  a  graph  of  &2  versus  r^/r.  for  various 
values  of  x. 
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Figure  8.  Bounding  Curve  for  Solutions  of  Equation  (60) 


F rom  Eq.  (62) 


a.  c>  <  a'  ■*  0  as  x  -»  0 
2s  2 


b.  a,  £  a'  -*  0  as  y  -*  x(u  -*  1) 

-3  4 


c.  For  any  u  =  —  >  1,  a^  is  maximized  for  x  =  1. 

i 

It  can  be  shown  that  a?  assumes  its  limiting  values  a_  =0  when  y  =  x. 

S  b8 

or  x  =  0. 

From  Eq.  (62)  it  follows  that,  when  x  =  1, 


cos 


ai  =  £* -  t  jm 

2  si  U3(l  +  u)  u  4  u 


iin  a'  =  J 


1  “  / 
1  -  cos  a' 


(63) 


Comparison  of  Eqs.  (33)  and  \ 6 3)  shows  that  a^  is  the  same  function  of  u  at 
x  =  1 ,  as  is  of  x.  It  thus  follows  from  Eq.  (34)  (also  see  Eq.  (44))  that 

a2g  <  <  K  =  5.  30°  (64) 

A  bound  can  be  placed  on  the  angle  eij  ,  which  is  similar  to  the  above  bound 

on  a ..  .  From  Eqs.  (56)  and  (58),  a,  <  ai  £  a.  ,  where  ai  is  defined  by 

2s  -i  \  >  ii  is  i  jni»  1  / 

the  following  equation  (see  Eq.  (61)): 

G(a')  =  maxK(a^)  =  H^H^. 
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or 


cos  a.j 


gpn  +  x) 

U2(u  +  X)2 


•i 


- 


x)(u2  + 


XU 


fin  =  >7 1  - 


2  , 
cos  a. 


(65) 


Figure  9  is  a  graph  of  versus  x  for  various  values  of  u  =  r^/r^. 
From  Eq.  (65)  and  the  previous  analysis,  it  is  clear  that: 


a. 

a.  £  a'  <  a 
Is  1 

b. 

Q1  =  alm  =  • 

c. 

aj  -*0  as  u- 

d. 

cij  =  0  when 

e. 

cij  =  arc  sin 

l 

u 


r. 

i 

rf 


Figures  10  through  16  were  prepared  by  Mr.  Jerome  Baker  using  the  results  of 
a  computer  program  based  on  the  technique  developed  in  this  report.  Figures 
10  through  12  show  the  minimum  bi -elliptic  velocity  increment  for  various 
values  of  rt  >  r^,  i.  e.  ,  for  intermediate  and  outer  bi-elliptic  transfers  (see 
Eq.(l)). 

The  plane  change  angles  0,  corresponding  to  Figures  10  through  12,  are  0*,  10°, 
20*,  and  30*.  respectively.  There  is  evidently  no  minimization  involved  for 
the  0=0®  case  which  is  treated  in  Reference  1. 


An  alternative  three-dimensional  transfer  procedure  is  the  modified  Hchmann 
transfer  described  in  Reference  2.  The  velocity  requirement  for  this  transfer 
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Figure  9.  Eounding  Curve  for  Solutionj  of  Equation  (59) 
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Figure  10.  Minimum  Bi- elliptic 


Figure  11.  Minimum  Bi-elliptic  Velocity  Requirements  When  0 
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Figure  12.  Minimum  Bi-elliptic  Velocity  Requirements  When  0  =  20 
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Figure  13.  Minimum  Bi- elliptic  v»lnr 
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Figure  15.  Division  of  Bi-elliptic  Plane  Change  for  Minimum  Velocity 
When  0  =  20*  =  a.  +  a,  +  a. 
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Figure  16.  Division  of  Bi -elliptic  Plane  Change  for  Minimum  Velocity 
When  0  =  30*  =  a,  +  a,  +  a 


mode  is  included  in  the  figures  for  comparison  purposes.  When  6=0*,  the 
modified  Hohmann  transfer  is  equivalent  to  the  simple  Hohmann  transfer 
(rt/r.  =  1). 

Figures  14  through  16  show  the  values  of  =  aj  8  and  which  yield 

the  minimum  values  of  AV^,  when  0  is  respectively  equal  to  10*,  20*,  and  30*. 
As  in  the  previous  graphs,  only  the  outer  and  intermediate  bi-elliptic  trans¬ 
fers  are  considered. 
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transfer  ellipse  to  a  given  transfer  altitude  hf  (radius  r^)  is  establis  ed.  When 
the  vehicle  reaches  hf,  a  second  impulsive  velocitv  change  AV 2  simultaneously 
changes  the  plane  by  amount  a 3  and  initiates  a  transfer  ellipse  to  the  altitude  hf 
(radius  r£)  of  tne  target  orbit.  A  last  impulse  AV3  changes  the  plane  by  amount 
Q3  and  circularizes  the  orbic  at  altitude  hf,  placing  the  vehicle  in  the  final 
(target)  circular  orbit. 

Studies  were  made  of  the  choice  of  plane  change  angles  a^,  and  a,, 
which  minimizes  AV  j  =  AVj  +  AV^  +  AV3  for  given  values  of  hf,  hf,  hf  ana 
total  plane  change  angle  0  =  a  j  +  +  Q3 .  Several  limiting  relations  were  ob¬ 

tained  lor  dj,  a^,  and  a^;  chey  are  dependent  on  either  the  ratio  rt/rj  alone,  or 
the  ratios  rt/r^  and  rt/rf,  and  are  independent  of  0. 
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